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The time dependence of quantum evanescent waves generated by a point source with an infinite 
or a limited frequency band is analyzed. The evanescent wave is characterized by a forerunner 
(transient) related to the precise way the source is switched on. It is followed by an asymptotic, 
monochromatic wave which at long times reveals the oscillation frequency of the source. For a source 
with a sharp onset the forerunner is exponentially larger than the monochromatic solution and a 
transition from the transient regime to the asymtotic regime occurs only at asymptotically large 
times. In this case, the traversal time for tunneling plays already a role only in the transient regime. 
To enhance the monochromatic solution compared to the forerunner we investigate (a) frequency 
band limited sources and (b) the short time Fourier analysis (the spectrogram) corresponding to a 
detector which is frequency band limited. Neither of these two methods leads to a precise determi- 
nation of the traversal time. However, if they are limited to determine the traversal time only with a 
precision of the traversal time itself both methods are successful: In this case the transient behavior 
^ ■ of the evanescent waves is at a time of the order of the traversal time followed by a monochromatic 

wave which reveals the frequency of the source. 
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I. INTRODUCTION 



> 
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' In order to summarize essential aspects of the time dependence of wave phenomena a number of characteristic 
velocities or times have been defined. The phase velocity, uj/k, is the velocity of constant phase points in the 
stationary wave (assume fc > for the time being) 



The boundary conditions, the superposition principle and the dispersion relation w — w(k) between the frequency 
Lo and the wavenumber fc determine the time evolution of the waves in a given medium. When a group of waves is 
y formed by superposition of stationary waves around a particular w, it propagates with the group velocity duj/dk. In 
^ ' dispersive media (where w depends on k), the group velocity can be smaller (normal dispersion) or greater (anomalous 
^ ' dispersion) than the phase velocity. It was soon understood that these velocities could be both greater than c for the 
propagation of light, and Sommerfeld and Brillouin |^], studying the fields that result from an input step function 
ILjl ■ modulated signal in a single Lorentz resonance medium, introduced other useful velocities, such as the velocity of the 
very first wavefront (equal to c), or the signal velocity for the propagation of the main front of the wave. Both the 
very first front and the signal velocity describe thus the causal response of the system and are therefore of particular 



H ' interest 



The above description is however problematic for evanescent waves, characterized by imaginary wavenumbers 
instead of the real wavenumbers of propagating waves. Their time dependence has been investigated by theoretitians 
and experimentalists in recent years because of its peculiar behaviour. A striking phenomenon is that, when crossing 
an evanescent region, certain initial wave features (such as the peak of the incident amplitude) appear at the far 
side at anomalously large speeds, but clearly comparison of peaks of transmitted and incident wave packets does not 
describe a causual process. Many publications and a recent workshop have been devoted to discuss the implications 
§. 

The role played by the imaginary part of the group velocity dw / dk and the possible definition of a signal velocity 
in the evanescent case have been much discussed but not yet completely clarified. Assume that a source is placed at 
X = and emits with frequency wq from t = on. If ujq is above the cutoff freguency of the medium (the one that 
makes fc = 0) a somewhat distorted but recognizable front propagates with the velocity corresponding to loq. Within 
the framework of the Schrodinger equation, and using a set of dimensionless quantities where the cutoff frequency is 
1 (see the subsection below), the dispersion relation takes the form |^ 

= 1 + fc^ (2) 

and the signal propagation velocity for the main front is equal to the group velocity, Vp — {duj/dk)^^ = 2{ljo — 1)^/^. 
In other words, at some distance x form the source, the amplitude behaves, in first approximation, as 
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(3) 



where fcg = (llJo ^ 1)^^^ is the wavenumber related to ujo by the dispersion relation, and O is the Heaviside (step) 
function. In the evanescent case, wq < 1, a preliminary analysis by Stevens [Q, following the contour deformation 
techniques used by Brillouin and Sommerfeld suggested that a main front, moving now with velocity Vm = 2(1 — 
luq)^/^ ~ lia{dLu/dk)^g, and attenuated exponentially by exp{nQx) (where Kq = (1 — ujoY^^), could be also identified, 

ip{x, t) « e"''^°*e-'^«^e(t - xv.m) ■ (4) 

The result seemed to be supported by a different approximate analysis of Moretti based on the exact solution 0| , and 
by the fact that the time of arrival of the evanescent front, r — x/v^m had been found independently by Biittiker and 
Landauer |^,^ as a characteristic traversal time for tunnelling using rather different criteria (semiclassical arguments, 
the rotation of the electron spin in a weak magnetic field, and the transition from adiabatic to sudden regimes in an 
oscillating potential barrier). 

In the treatment of Stevens, as well as in the original work by Sommerfeld and Brillouin, the contour for the integral 
defining the field evolution was deformed along the steepest descent path from the saddle point; and the main front 
was associated with a residue due to the crossing of a pole at ikq by the steepest descent path. But later, more 
accurate studies of the punctual source problem or other boundary conditions showed that the contribution from the 
saddle point (due to to frequency components above or at the frequency cutoff created by the sharp onset of the source 
emission), and possibly from other critical points (e.g. resonance poles when a square barrier is located in front of the 
source [^) were generally dominant at r, so that no sign of the a>o~front (^ can in fact be seen in the total wave 
density at that time fq-p^. Similarly, corrections to the original work by Sommerfeld and Brillouin have been also 
worked out for electromagnetic pulse propagation [ p^ . In spite of these clarifying works, several important aspects 
have remained obscure or not investigated, such as the actual time scale for the attainment of the stationary regime, 
the characterization of the transients, and the role (if any) played by t in the time dependence of the quantum wave. 

Recently, one of the authors in collaboration with H. Thomas [Q, reconsidered the problem of Sommerfeld and 
Brillouin, and provided a detailed discussion of the forerunners and the signal sent out by a source which has a sharp 
onset in time. These authors also pointed out that the forerunner, generated by switching on the source, is associated 
with a time-dependent wide band spectrum whereas the signal, the wo^front, carries the oscillation frequency of the 
source into the evanescent medium. The signal is called a " monochromatic front" . In contrast to previous work, which 
tried to find a front simply by analyzing the amplitude of the waves, these authors emphasized the frequency content 
of the forerunner and the signal. In the evanescent case, the amplitude of the monochromatic front is exponentially 
small compared to the forerunner, and in agreement with the works mentioned above, it cannot be detected using a 
simple criterion based on the magnitude of the wave. Two approaches were proposed to enhance the monochromatic 
fronts compared to the forerunners. First, the dominance of the forerunners might arise due to the fact that high 
frequencies are transmitted in the propagating energy range. This can be avoided if the source is frequency limited 
such that all frequencies of the source are within the evanescent case. (Technically this means that a frequency 
window is chosen to avoid the effect of the saddle point contribution.) A second option is not to limit the source 
but to frequency limit the detection. We can chose a detector that is tuned to the frequency of the source and that 
responds when the monochromatic front arrives. 

The aim of this work is to characterize the time dependence of Schrodinger evanescent waves generated by a point 
source. Wc identify several wave features, in particular the arrival of the first main peak and the transition from a 
forerunner dominated behavior to an asymptotic behavior dominated by the monochromatic front. We also investigate 
in some detail the proposals made in Ref. |14| to enhance the monochromatic front and consider both frequency limited 
sources and a frequency-time analysis of the wave at a fixed position. This leads to the investigation of the spectrogram 
of the wave generated by the source. 

For a source with a sharp onset, we find that the traversal time r plays a basic and unexpected role in the transient 
regime. For strongly attenuating conditions (in the WKB-limit) the traversal time governs the appearance of the 
first main peak of the forerunner. In contrast, the transition from the forerunner to an asymptotic regime which is 
dominated by the monochromatic signal of the source is given by an exponentially long time. If the source is frequency 
band limited such that it switches on gradually but still fast compared to the traversal time, the situation remains 
much the same as for the sharp source, except that now the transition fom the transient regime to the stationary 
regime occurs much faster, but still on an exponentially long time-scale. The situation changes if we permit the source 
to be switched on on a time scale comparable to or larger than the traversal time for tunneling. Clearly, in this case 
a precise definition of the traversal time is not possible. But for such a source the transtion from the transient regime 
to the asymptotic regime is now determined by the traversal time. Much the same picture emerges if we limit instead 
of the source the detector. As long as the frequency window of the detector is made sharp enough to determine 
the traversal time with accuracy, the detector response is dominated by the uppermost frequencies. In contrast if 
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the frequency window of the detector is made so narrow that the possible uncertainty in the determination of the 
traversal time is of the order of the traversal time itself, the detector sees a crossover from the transient regime to the 
monochromatic asymtotic regime at a time determined by the traversal time. 

Possibly, the fact that we can not determine the traversal time with an accuracy better than the traversal time 
itself tells us something fundamental about the tunneling time problem and is not a property of the two particular 
methods investigated here. 

A. Dimensionless quantities and notation 

The (dimensional) time dependent Schrodinger equation for a particle of mass m moving in a constant potential 
V{X) — V is given by 

The number of variables and parameters may be reduced by introducing dimensionless quantities for position, time 
and wave amplitude, 

h 

TV 

This allows to write the corresponding dimensionless Schrodinger equation, 

■ dtp <9V , , 

Other useful dimensionless variables related to the dimensional energy E (or frequency W = E/h), and wavenumber, 
K — [2m{W — V/h)/h]^/'^, are respectively 

Lu = E/V (10) 
The reader may check that the dimensional dispersion relation 

takes for dimensionless quantities the simple form given in (^). 

II. SOURCE WITH A SHARP ONSET (INFINITE FREQUENCY BAND) 

In this section we shall investigate the time dependent wave function corresponding to the "boundary condition" 

il;{x = 0,t) = e-"^°'e{t), (13) 
which may also be given by the corresponding Fourier transform 

ip(^x ^0,uj) ^ -. (14) 

The superposition 
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■ poo ^ikx — iujt 

i,{x,t) = — dio — - (15) 



satisfies the Schrodinger equation as well as the boundary condition (|13[). Along the integration path, k is positive 
for Lu > 0, and purely imaginary (with positive imaginary part) for w < 0. This corresponds to outgoing waves 
from the source. In Ref. p4| the integration was carried out in the complex frequency plane since this permitted a 
close comparison between the calculation for the Schrodinger equation and for relativistic field equations. If only the 
Schrodinger equation is of interest the complex fc-plane is most advantageous. In fact this permits to express the 
integral in terms of known functions. In the complex fc-plane we have, 



^{x,t) = ~— dk2k , (16) 



2-1X1 



dk 

+ 
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fc + ikq k — iKQ 



where the contour goes from —oo to oo passing above the pole at ikq- The two terms in (|lj) lead to integrals 
with the form discussed in the Appendix A. The contour can be deformed along the steepest descent path from the 
saddle at kg = x/2t, the straight hue 

ki = ~kR + x/2t, (18) 

(kn and fc/ are the real and imaginary parts of fc.) plus a small circle around the pole at iKo after it has been crossed 
by the steepest descent path, for fixed x, at the critical time 

r^:^- (19) 
2kq 



This procedure allows to recognize two w- functions |15|Jl6| | (see the Appendixes A and B), one for each integral, 

i;{x, t) = ie-*+^'^^ [w[-u',) + w{-u'^)] . (20) 

Here, 

-0 = ^i^/^-o {-^ - \) (21) 

It is clear from the exact result (^0[^l]), that r is an important parameter that appears naturally in the w- function 
arguments, and determines with the global properties of the solution. Its detailed role will be discussed in the 
following sections . 

Eq. ( pO| ) is in agreement with a previous expression by Moretti, derived using different contour deformations and 
notation. Our analysis of this exact result will be however quite different and more detailed. 



A. Approximations 

Often an exact expression is not very informative by itself, and the approximations make its essential content 
manifest in certain limits. The simplest approximation for 'ip{x,t) for times before t is to retain the dominant 
contribution of the saddle [by putting fc = fc^ in the denominators of ( p!7[ ) and integrating along the steepest descent 
path]. It is useful to write this in different ways, 

g-zt+ifc^ 1 _ e-'*+'^'V(2t/^)i/^ 

- 2z7ri/2 [y^^VJ- (z-l)«o(T2+i2) 

27r V Vlst fls-Qo ' 

where 
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17, = fcf = — , no ^-4 (23) 

(These are particular values of the frequency variable = w — 1, defined with respect to the frequency level of the 
potential.) The average local instantaneous frequency for this saddle contribution [Q is equal to the frequency of the 
saddle point, 

ujs^l + fls. (24) 

(Note the different sign of fl^ in the exponent of Eq. (p^). Qg depends on and a time derivative has to be taken 
to obtain (||), see Q.) 

After the crossing of the pole ikq by the steepest descent path at < = r the residue 

V'o(a;,i) = e-^'^°*e-'=°^e(i-T). (25) 

has to be added to (p^, 

t^ix,t)^iP,{x,t)+Mx,t). (26) 



The solution given by Eq. (^5|) describes a monochromatic front which carries the signal into the evanescent medium. 
The conditions of validity of this approximation can be determined by examining the asymptotic series of the w{z) 
functions in ( po|) for large \z\, see the Appendix B. In fact ( p6| ) is obtained from the dominant terms of these expansions. 
The modulus of the two w— function arguments in (^) is given by 

Mit)^\~u'o\^\-~u'^\^^{1^+ry/\ (27) 

This quantity may be large in different circumstances. It goes to oo when kq, t, or x go to oo, and also when t —^ 0. 
M{t) is minimum at the crossing time t = t, where it takes the value M(r) — {xkoY^^- For xk < 1 the pole lies 
within the range of the saddle Gaussian and cannot be treated separately. 

Assumming that M is large the phase of z determines the appropriate asymptotic expression of w{z). For Im{z) > 0, 
w(z) ~ i/{-K^/'^z)^ whereas for Im{z) < 0, w{z) ~ i/{T:^^'^z) + , see (B5) and (B6). arg(— Uq) goes from 7r/4 to 



Stt/A when t goes from t = to oo, and arg(— Mq) goes from 7r/4, at t — 0, clockwise to — 7r/4 when t oo, and 
crosses the real axis at t = t. Thus the application of the previous asymptotic formulae lead exactly to Eq. (p6|). As 
discussed in the Appendix B one should not be mislead by the seeming front in (p6|). It is possible to obtain a smooth 
approximate expression around t = r by adding a correction to ^ps that takes into account the region near the pole 
This corrected expression however is asymptotically equivalent to ( |2^ ) because at r, and within the conditions 
that make the saddle approximation valid, the contribution of the pole is negligible. To see this more precisely let us 
examine the ratio between the modulus of the two contributions, 

\lps\ X 

Its value at r is an exponentially small quantity, 

R{t) = e-'""'{2TTKoxy/'^ . (29) 

Note also that R has a minimum at r/S^/'^ < r and tends to oo as t — s- oo. As a function of x, the minimum of R 
decreases up to xkq = 1/2, and then grows again monotonously. 

In summary, for the source with a sharp onset described here, the monochromatic front is not visible when the 
approximation (|2^) remains valid around t — t. However two very important observable features of the wave can be 
extracted easily from (p6|). The first one is the arrival of the transient front, characterized by its maximum density 
at tf = It is important to emphasize that this time is of the order of r, but the wave front that arrives does 

not oscillate with the pole frequency ujq, but with the saddle point frequency uJs- Figure ^ shows this transient front 
for three positions. In this and similar figures the densities are exponentially amplified by e^**"^ to make possible the 
comparison among different values of x with the same scale. Thus, all amplified densities A = \ip\'^e'^'^°^ tend to one 
in the asymptotic large t regime, and the corresponding logarithm to zero. Also shown is the approximation due to 
ips, althout it is only distinguishable from the exact result for x = 7. Note that the amplified density j-^spe^""^ is 
simply R~^ . In Figure ^ the instantaneous average frequency w is represented for the smallest x value of Figure 
X = 7. At small t and around the transient front, Hj — 
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The second observable feature that we can extract from ( pq ) is the time scale for the attainment of the stationary 
regime, or equivalently, the duration ttr of the transient regime dominated by the saddle before the pole dominates. 
ttr can be identified formally as the time where the saddle and pole contributions are equal, R — 1. Because of ( |29| ) 
we shall assume t « ttr to obtain the explicit result 



2 /3 



2K27rl/2 



(30) 



Using kqx >> 1, the velocity for the motion of the the space "point" where the transition from transient to stationary 
behaviour (or from saddle to pole) takes place is given by 

Vtr ~ ■ (31) 

Contrary to the motion of the front maximum, this point does not move with constant velocity. The transition may 
be observed in various ways. In Fig. || the logarithms of the exact densities and of their components are represented. 
The crossing point where R = 1 may be observed for a; = 7 in the change of behaviour of the total wave density. The 
oscillation around ttr is due to the interference between -00 and ^s, and has the characteristic period Tq = 27r/|r2o|- 
In Figure I the crossover between the regime dominated by the saddle frequency and the one dominated by loq is also 
easily noticeable. 

Finally, when xhq is small 1), the saddle approximation describes correctly the very short time initial growth, 
but fails around r because the pole is within the width of the Gaussian centered at the saddle point. The pole cancels 
part of the Gaussian contribution so that the bump predicted by V's at r/S^/^ is not seen in this regime. Fig. ^ shows 
the time dependence of the density for two small values of x. r does not correspond to any sharply defined feature, 
but provides here a valid rough estimate of the attainment of the stationary regime. 



III. FREQUENCY-BAND LIMITED SOURCE 



H. Thomas and one of the authors suggested recently to avoid the dominace of the saddle point solution over 
the monochromatic front by limiting the frequency band of the source. Specifically, it was proposed to cut-off the 
w-amplitude (|l5). 



ip{w, X — Q) 



1 



— - [loo - Auj)) - e{uj -{uj- Auj))] 



(27r)i/2cj-cjo + i0 

This implies that the emission of the source is not sharply defined as in (M), see the Appendix C. Instead, 



ipix = 0,t) 



e{t) - -ImEii-iAujt) 



(32) 



(33) 



where Ei is the exponential integral defined in (C4). The frequency band width Ao; may be chosen so that the onset 
of the source is fast with respect to t and that all the frequencies in the source are in the evanescent region [p^ , 



27r , 

T 



(34) 



Later we will see that in fact it is necessairy to give up the condition that the source switches on fast compared to 
the traversal time. Now "negative times" are also required to describe the signal growth. For an arbitrary x, 



ipix,t) 



duj ■ 



-'i{ijjt—kx) 



271" Juo-Auj uj-ujQ + iO 



(35) 



In this case the integration technique used in the previous section does not provide an analytical solution. We shall 
manipulate the integral to facilitate the exact numerical evaluation and the discussion of approximations. Using 
— UJ — 1, Eq.(p5|) becomes 



ip{x,t) = 



it /.r2o+Aa) 



27r Jn„-Au, n-flo + iO 



(36) 
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The integrand of Eq. ( p6 
chosen according to Eq 



has a saddle point at fis = (a:/2t)^. But since JIq is now negative, a frequency interval 
3J), now excludes this saddle. It is convenient to introduce the variable y = U,/VLs so that 



the valley-hill structure of the exponent remains constant, 

jg-»* ry+ giAsign(t)(2i,i/^-y) 

^{x,t) = ^ dy — — . (37) 

27r Jy_ y-yo + iO 

Here, 

yo = i^o/^s , y± = 2/0 ± Au/ils , A = |risi| . (38) 

In order to perform contour deformations in the complex y plane it is also necessary to specify that the branch cut for 
j/^/^ is set along the positive real axis, and that the saddle is at y = 1 for t > (just above the cut) but at ?/ = e^*'^ 
for t < (just below). The steepest descent and ascent paths from it are given, in the first Rieman sheet of y, by the 
sections of the parabolas 

= sign(i)^(l - 2/fl)2 (descent) (39) 
2// = -sign(t)i(l - yj^)^ (ascent) (40) 



starting from the saddle point. Irrespective of the sign of t the original contour in ( pT[ ) is entirely within one of the 
valleys of the saddle, so that there is no need to take this critical point into account in the contour deformation. The 
most efficient contour deformation consists on following the steepest descent path from the lower integration extreme 
2/_ downwards to infinity and coming up to the upper integration extreme y+ following its steepest descent path. If 
t > this contour encloses the pole at loq and the corresponding residue has to be included. The extreme points 
become the critical points of the integral, apart from the pole at uq when t > 0, 

4'ix,t)=D^-D++^b'o, (41) 

where D± (generically D^) are the integrals from y± to oo along the corresponding steepest descent paths, SDP{z), 



and 



D, = ^ dy , (42) 

27r JsDP(z) y-yo 



i^(,(x,t) =e-*"«*e-"«"e(i). (43) 



Unlike ( pSf ) the residue is now present at all positive times. Using the generic notation z — y± for any of the two 
extreme points the steepest descent paths from them and their slopes are given by 

yi = isign(i)[y4 _ Aylil + z) + 2y^z(4 + 3z) - Az^yR{l + z) + z^]'/^ (44) 



dyi ^ yl ~ iylil + z) + ynzjA + 3z) - z^(l + z) 
dyn 2yi 
dyi 



(45) 



= sign(i)(~z)i/2. (46) 

oyB^ z 

In the numerical evaluations the value of y along the path may be obtained by solving for each step 

dy ^ dyi{dyR/dyi + i) , (47) 

with the initial condition y — z. 
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A. Approximations 



The leading term in A^^ of ( p^ ) is obtained by integrating from z along the ray with the direction of the steepest 
descent path at z and taking all functions in the integrand, except the exponential, out of the integral, with their 
values at z ]l7t , 

^ Dl ^ —- \ -~u..~...^.0. (48) 

2nX{z ~ yo)\z - 1| 

where lUz = w± = ± Aw, Kz = k± — {1 — ^jJ±Y^'^ , and 6z is the angle of the steepest descent path at z, 

0z=Tr + arctan[sign(i)(-z)^/2] (49) 
(the branch of the arctan function between —tt/2 and 7r/2 is taken). The modulus of these contributions is given by 

'^^'= 2.A.y+4)^/- ' ^''^ 

where 

are the critical times when the steepest descent path from the saddle crosses the extreme points of integration y±. 
Regarding the possibility of making the monochromatic front visible, the frequency band limited source elliminates 
the saddle dominated transient effects, but they are substituted by transients associated with a new critical point: 
the upper extreme of integration. The actual time for transition to the stationary regime is still larger than r under 
semiclassical conditions, but much smaller than for the sharp onset case. A difference with the infinite frequency 
band case is that the maximum of \D±\ is at i = 0. This is quite remarkable. A particular feature of the wave, its 
maximum, appears instantly with zero delay at arbitrarily large distances. One should keep in mind though that the 
build up of the initial signal has required previously an infinite time, from t — ~oo to t — 0. 

The transient regime previous to the dominance of the pole must be also different from the infinite band case 
since now there is no saddle point contribution. The transient is here dominated by the integral D+ because of the 
exponential dependencies and by the frequency of the upper integration limit An analysis similar to the one 
performed in the infinite frequency band case can be now performed. Let us compare the contributions from the pole 
and by dividing their moduli, 

^, ^ ^ ^ e(—o+^^)-[2nALoit^ + tl)'/^] . (52) 
\D+\ 

From the condition i?' = 1 we obtain the new critical time 



„2x(ko-k+) \ „£c(ko-«;+) 

and the critical velocity 



~ \ 4^2 A^2 4) ^ 27rAa; ' ^^^^ 



(54) 



They have to be compared with the corresponding quantities for the sharp onset case, (|3^) and ( pT| ) respectively. The 
critical time for the transition depends exponentially on x in both cases, but it arrives exponentially earlier for the 
frequency band limited case. In Figure ^ the logarithm of the amplified density and of the approximation provided 
by |i?1|2 is represented for two values of x. For a; = 50 the transition time ttr is relatively small and can be seen 
in the time interval shown (The oscillation period of the interferences around ttr is now 27r/Acj). For x = 135, Ur 
is much larger and cannot be seen in the scale chosen; note the well defined peak at t = 0. The slight disagreement 
between the exact result and the approximation provided by is because, for very small times, the slope of the 
steepest descent path from y+ is very small and the path passes close to the pole. The pole perturbs the integral in 
this manner for times such that A^^ < Aui/ils, i.e., t < 27r/Aw. 
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Prom Eq. ( p^ ) we see that the dominance of the monochromatic front at t and fixed x requires that the magnitude 
of (— Ko + K+)x ft! 1 or smaher. This can be achieved in two ways: (a) We can consider a frequency wq, and a 
frequency interval such that — kq and k+ become small. This means that both 1 — ti^o and Aw should be close to 
zero, but to avoid propagating waves Aijj < (1 — wp); (b) Alternatively, we can allow ldq to take any other value below 
1, but making the the frequency interval Aa; small. In both cases a small frequency window is needed and we can 
expand with respect to Aw so that (— + k.+ )x = — Awr. Thus we see that the attempt to measure the traversal 
time accurately has to be abandoned. The lower limit of Eq. (|3j) is now violated. Figure ^ shows an example where 
the first inequality in (|3^ ) is not obeyed. The traversal time (rather than ttr^ which in this case is an imaginary 
number) now marks the transition to the asymptotic region even though the build up time is so large that the arrival 
of the monochromatic front cannot be defined with a precission better than r itself. In this respect. Figure ^ is very 
illustrative. It shows how the transition of the average frequency characterizing the transient to the frequency of the 
monochromatic solution occurs around the traversal time r, but the time interval required for the transition is clearly 
larger than r. 



IV. TIME-FREQUENCY ANALYSIS OF THE WAVE FUNCTION 

Instead of limiting the source in frequency we can limit the detection of the field to a range of frequencies of 
interest [14]. The wave amplitude at fixed position becomes a function of time, i){t), that can be Fourier analyzed 
to provide its frequency representation, ip(uj). However, neither of the corresponding densities, [^'(OP o^' 
tells when a particular frequency arrives or decays, nor what is the relative importance, at a given time, of different 
frequency components. This type of information is provided by joint time-frequency representations p3|l . There 
are many possible ways to carry out a time-frequency analysis. A simple one, and surely the most common, is the 
"spectrogram" based on the short time Fourier transform (stFt). We are interested in the Fourier spectrum that 
would be measured at the observation point x if the wave is observed during a time interval of duration T . If a precise 
determinantion of the traversal time is attempted, the extend of the time interval T is chosen to be short compared to 
the traversal time t. On the other hand, if T is too short then the frequency resolution will be poor. A compromise 
to obtain a good resolution in time and frequency is for T to be bound above and below [14], 

?^£.i«i«E^.|l. (55, 

X T 1 111 In 

This is equivalently expressed in terms of inequalities for time scales 

r » T » To, (56) 

where Tg = 27r/|ilo| is the oscillation period corresponding to the frequency |r2ol- Combining the inequalities one 
finds K,[)X » 1, namely when these conditions are satisfied the saddle contribution to the wave function will be a 
very good approximation to the total wave up to an "exponentially long" time Ur » t (see section 2). 
The short time Fourier transform (stFt) of the field ip{x, t) is given by 

I rt+T/2 

F{uj;x,t) = -—— di'exp(iwt')^(a:,t')- (57) 

(27r)V-: 7t_-r/2 

Note that the short time Fourier spectrum depends, in addition to the frequency, parametrically on the observation 
point X and the time t. 

Consider now first the contribution of the monochromatic front to the stFt. We denote this stFt by Fp{LLj;x,t). For 
t <T~ T/2 we have Fp(w; x, t) = 0. For t > r + T/2 we find 

FA^;.,t) . -^e-»-e-(--o)- -^-[(— (58) 
(27r)V^ [lo - LUo) 

This amplitude peaks at the frequency of the source ujq and, for t > t + T/2, it is time-independent at this frequency. 
Away from this frequency the stFt of the monochromatic wave decays algebraically and oscillates sinusoidally. In the 
time interval r — T/2 < t < t + T/2 the stFt tracks the arrival of the monochromatic front, 

F>;.,t) = ^ ZrrZT, • 
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Consider next the saddle point contribution. Its frequency in a short time interval is determined by the expansion of 
^ls{t')t' away from t. This expansion gives ^ls{t')t' = 2^ls{t)t — il.s{t)t' + 0{{t' — i)^). Taking into account that the 
prefactors are slowly varying over the time interval of interest here (this is assured by the first inequality in (|55|)), we 
obtain, 



The stFt of the saddle point solution peaks at the frequency tUs = + 1. This frequency is very large at short 
times ("kinetic regime"), where ~ = a;^/(4t^) is dominated by frequencies above the potential. At long times 
("potential regime") LUg tends to the potential frequency 1. The transition time between these two regimes of different 
ujs behaviour may be estimated, by solving fls{t) = 1, as t — x/2 (it can be larger or smaller than r.) 

We see that the two Fourier transforms peak at well separated frequencies. Still the stPt of the pole has an 
exponentially small amplitude compared to that of the forerunner. Therefore, it is possible that the stFt of the saddle 
is still large at the frequency ojo where the stFt of the monochromatic front peaks. Thus we have to investigate the 
StFt of the saddle at the frequency luq. From Eq. (|60| ) we obtain. 



At a time t ~ t the amplitude of the saddle is still of order 1. Thus even at the peak frequency of the front its 
contribution to the spectrum is of order 1 and much larger than the exponentially small peak of the monochromatic 
front. 

Let us now investigate the properties of the transient in the frequency-time domain, and determine the role played 
by r, in particular at the frequency of the signal loq. The spectrogram is defined as the square modulus of the stFt, 
S(t,uj;x) — N\F(uj; x,t)\'^ , where TV is a normalization constant. (The notation of the argument of S is appropriate 
for a time-frequency analysis at fixed x.) For "normalizable" cases N is chosen so that / J dtduS — 1. In our case, 
this is not possible because the asymptotic stationary regime does not decay, but S provides anyway information on 
the relative importance of two time-frequency points, so N is chosen to be some convenient value, N = 7r^/2. 

For analyzing the transient we may neglect the monochromatic front and concentrate on the saddle point term. 
The spectrogram of the saddle point solution is denoted by Ss{t^uj;x) and is given by 

Ss{t,u;;x)=a{t) 62 



with an amplitude 



Note that a{t) is also proportional to the absolute square of the saddle point solution \tps{x,t)\'^ — (4/7r)Q(f). The 
amplitude a has a maximum as a function of time at the transient front peak, tj — S^^^^r, that moves with a speed 



\/3«™. (64) 



Thus Vf determines both the speed of the peak value of the saddle point solution |i/)s(a;,t)p and the speed of the 
absolute maximum of the spectrogram. The peak value of the spectrogram for a given time is at the frequency 
= fis -|- 1, see an example in Figure |^. Consider next the spectrogram at the frequency of the source, u = ujq. For 
fixed X it is bounded by the envelop function 

and has local maxima close to the times at which the sin function is ±1. I3{t) grows as for short times and decays 
as t~'^ for large times. Its maximum arrives, for fixed x, at -^/S/Sr. Equivalently, this maximum moves with a speed 

Venv = \/3/5Vm,- (66) 

This reveals that, in 5(^,^0; x), r is not the time where the stationary regime begins (it is necessary to wait an 
exponentially long time to attain that regime), but the basic time scale for the arrival of the main part of the 
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transient. The detailed oscillatory pattern of this main part will depend on the value of T /Tfj, see the Figures || and 
^. For to = loq the sin function vanishes at times 

" (nTo/T-l) 

(In the uj ~ t plane the spectrogram has a maximum which follows the line uj ~ + 1 and vanishes along the lines 
Lo = Vlg + \ + 2TTn/T. In between these zero-lines the spectrogram exhibits local maxima.) The n for the closest zero 
to r is n w 2T/Tq whereas the largest i„ (associated with the minimum n so that (|6^) has real solution) corresponds 
to n w T/Tq. Hence there are approximately T/Tq zeros between r and the last oscillation. For T/Tq < 1 there is a 
single major bump close to r, and for T/To >> 1 there are many local maxima that "sample" the form of the envelop 
function (|65|). 



A. Other time-frequency distributions 



The time-frequency characterization of the wave function at fixed position depends on the quasi-distribution chosen. 
There is nothing wrong with this non-uniqueness as long as the quasi-distribution chosen is specified. In choosing the 
distribution one may consider several factors. Ideally the distribution should be easy to calculate, and should not be 
too noisy (e.g. with wild oscillations). It may also occur that one of the distributions is naturally adapted to the way 
the detection experiment is performed. 

Let us first discuss a different type of spectrogram. The roles of time and frequency can be inverted so that instead 
of limiting the Fourier transform of ^plt) with a "square window" 

h{t) = Q{t-T/2)Q{-t-T/2), (68) 

a "short frequency Fourier transform" (sfFt) may be similarly defined for by limiting the frequency integral with 
a window, 

g(w) = e(w - Acj)e(-w - Aw). (69) 

The corresponding spectrogram, 5", is obtained as the square modulus of the sfFt. The two spectrograms, S and 5', 
are however not equal because the window functions h and g are not Fourier transforms of each other. In particular, 
note that the spectrogram S' may be equivalently obtained from a short frequency Fourier transform analysis for an 
initially sharp onset signal, or as the density that results from a source with a smooth onset (The density of ip{x,t) 
in Eq. ( ^5| ) is proportional, up to a trivial normalization constant, to S' {t,uJo; x).) As a consequence, S'{t,ujQ) peaks 
at t — after an infinite time growth whereas S{t, luq) is strictly zero at t < ~T/2 and its main transient part 
peaks around r. The qualitative behaviour is different because of the different window functions, even though the 
inequalities used here or in section III, Eqs. ( [ss] ) and (^4|), are in fact identical if one identifies T — 2tt/ Auj. (Compare 
Figure ^ and the upper curve of Figure ^ where the values of r, T and Tq are equal.) These inequalities imply in both 
cases that the pole contribution can be entirely neglected at or around r, and that it will only be of importance at a 
a time that depends exponentially on kqx. 

We have also considered the time-frequency Wigner function. This representation separates much more clearly the 
saddle and pole contributions (it also associates to the later an initial time t = rather than t = t), but presents the 
inconvenience of a very rapidly oscillating pattern and strong interference terms. 



V. SUMMARY AND DISCUSSION 



The time dependence of evanescent Schrodinger waves created by a point source has been investigated, combining 
analytical or numerical exact results and approximate expressions. The background of the former allow to test and 
contrast the validity of the simplified description of the later. We have also performed a time-frequency analysis of 
the transients which preceed the stationary regime by means of spectrograms. 

An important aspect of the work is the elucidation of the role played by the different parameters, in particular 
by the "traversal" time r. It is a basic parameter that determines the global shape of the wave but at this time, in 
semiclassical conditions kqx > 1, the monochromatic front (with the main frequency of the source) is not observed in 
the total wave because of the dominance of the saddle point contribution. However, for the source with a sharp onset 
we have identified a transient front whose maximum peak arises at a time The dominant frequency of this 

front does not correspond to the main frequency of the source, wq, but to the mean saddle point frequency (above or at 
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the cutoff frequency). The time ttr > t that determines the duration of the transient, or the passage to the asymptotic 
regime, has also been identified. Frequency band hmited sources accelerate the crossover to the stationary regime but 
there is still a transient dominated by the upper frequency contained in the band. The spectrogram S{t,uj;x) has 
been also investigated to show the variation with time of the frequency components. In this representation the main 
contribution at the frequency ujq arrives at x around the traversal time r. This main contribution is a transient due to 
the saddle point. At an exponentially long time the pole term will eventually take over and remain as the dominant 
contribution in the stationary regime. 

In fact it is possible to see the mochromatic wave front arriving around the traversal time r, but only when the 
semiclassical conditions are abandoned, and with an accuracy which is never better than the traversal time itself. 
This occurs if the point source starts the emission abruptly, but also when it is switched on gradually (limiting the 
frequency band), or for a limited frequency band detector, whose response is modelled here with a spectrogram. The 
coincidence of all these cases strongly suggests that the limitation of the accuracy to measure r is in fact a general 
property. 
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APPENDIX A: INTEGRALS 

In this appendix we shall solve the integrals of the form 

1= / dk— ; a>0, (Al) 

where r_|_ goes from —oo to oo passing above the pole at fc = fco- The saddle point of the exponent is at A; = —b/2a and 
the steepest descent path is the straight line kj = —{kii + b/2a). The original contour is at the border between hill and 
valley and can be deformed into this path, taking into account the residue of the pole when Im(fco) < — [Re(fco) + &/2a]. 
By completing the square and introducing the new variable 

M = u(fc) = ii^ai/2(fc + &/2a), (A2) 

the integral takes the form 



I^e'^ { / du - 2i7re-"oe[Im(uo)] } (A3) 

[J-ooU-uo J 

where uq = u{k = kg). Note tha t th e m— varia ble has its origin at the saddle point and its real axis corresponds to 



the steepest descent line. Using (B2) and (B3) it can be finally written as 

- 

X = -iiTe'~ w{-uq) . (A4) 

APPENDIX B: PROPERTIES OF w{z) 

The iti— function p^ , p^ is an entire function defined in terms of the complementary error function as 

2 

w(z) = erfc(-zz) . (Bl) 
w{z) is frequently recognized by its integral expression 

1 /■ e""'' 

wiz) = — / du (B2) 

in Jy_ u — z 
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where r_ goes from — cx) to oo passing below the pole at z. For \ra.z > this corresponds to an integral along the real 
axis. For Imz < the contribution of the residue has to be added, and for Imz = the integral becomes the principal 
part contribution along the real axis plus half the residue. From (B2) two important properties are deduced, 



and 



■w{—z) = 2e ^ — w{z) 
wiz*) ^ [wi^z)]*. 



(B3) 
(B4) 



To obtain an asymptotic series as z — *■ cxd for Imz > one may expand (u — z) ^ around the origin (the radius of 
convergence is the distance from the origin to the pole, |z|) and integrate term by term. This provides 



w{z) 



^ 1-3. 

rn— 1 



(2to- 1) 



{2z^y 



Imz > 



which is a uniform expansion in the sector Imz > 0. For the sector Imz < (B3) gives 

'w{z) ~ 



^ 1 • 3 • ... • (2to- 1) 

m— 1 



(2z2)" 



2e- 



Imz < . 



(B5) 



(B6) 



If z is in one of the bisectors then — z^ is purely imaginary and the exponential becomes domin ant . But righ t at 
the crossing of the real axis, Imz = 0, the exponential term is of order o(z~"), (all n), so that (^5| ) and (B6) are 
asymptotically equivalent as |z| oo. One may insist however in removing the discontinuity of the series expansion 
(as a function of z) around the real axis. This can be done by adding a correction to the saddle contribution that 
takes into account the region close to the pole in the line integral, see also the discussion in Rcf. |14 . Note that in 
(B5) only the region around the origin contributes. Consider now also the region around the pole. 



-du ' 



Zr + A 



ZR-A 



du 



In 



zr - A - z 



(B7) 



For large |z| and in the proximity of the real axis this may be approximated by e ^ Z7rsign(z) which exactly cancels 
the discontinuity between (B5) and (|Bq). 



APPENDIX C: TIME DEPENDENCE OF FREQUENCY BAND LIMITED SOURCE 

The time dependence of the frequency band limited amplitude at x = is obtained from ( ^^ by Fourier transform 

fuiQ+Aui ^^ — iiAjt 2 



4^{x — 0,t) = — I duj ■ 



27r Jujo-Auj uj~ujQ + iO 



—e 



-iujQt 



^ pi^o+Auj -—iuji 



Using the frequency uj' = (uj — (x'o)i, 



UJQ + Au} 



duj ■ 



UQ-Auj w - Wo 



tAuj -iu 

du;' — 



(CI) 



(C2) 



and the principal part integral can be expressed by contour deformation in terms of combinations of exponential 
integrals |p5||, 



dY ■ 



Y 



El {~ib) — El (ia) — ztt , a, 5 > , 



where 



Eiiz) 



dY- 



Y ' 



argz < TT 



(the contour does not cross the negative real axis). Finally, using (|C3|), (CI) can be written, for all t, 



(C3) 



(C4) 



(C5) 
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APPENDIX D: RELATIVISTIC CASE 



Defining the dimensionless parameter c by combining the dimensional velocity of light C, the mass and the potential 
constant, as 



c = C(2m/y)i/2, 

the dispersion relation for a Klein-Gordon wave equation takes the form 



(Dl) 



(D2) 



For evanescent conditions k is purely imaginary, k — in. In particular for lo = luq, or {ft — fio), kq cannot take 
arbitrarily large values, < kq < c/2, contrast this with the non-relativistic case where there is no upper limit. 
For a source with a sharp onset the wave function is Mm 



fp(x,t) = — e" 

ZTT 



1 



which is zero for t < x/c. There are now two saddle points at 



±f7, = ±- 



29 



(D3) 



(D4) 
(D5) 



and two branch cuts. The integration contour may be deformed along the two steepest descent paths. The pole is 
crossed at |14l 



X 

2ko' 



(D6) 



Note the lower bound x/c < t when 17 — > (or kq — > c/2). Consequently the traversal time is strictly limited by the 
velocity of light Q . 

We shall only discuss the contribution from the saddle at +ris corresponding to the excitation of particles , 



1/2 



e{ct-x). 



(D7) 



Line arizing the square root of the exponent to evaluate the short time Fourier transform of the saddle wave function 
(D7), one obtains 



Ss(t,Lu;x) 



2x' 



nl sin2[(f7 - f7o)r/2] 



In this case the maximum of the envelop as a function of t is given by 



[t > T/2 + x/c). 



T 
3172 



1 + 3 



"0 
c2/2 



-I 1/2 



(D8) 



(D9) 



Because of the dispersion relation Qa < c^/2 in the evanescent case, it follows that te,i is bounded by and 
(4/3i/2)r. 



[1] L. Brillouin, in Wave propagation and Group Velocity, (New York, Academic Press, 1960). 

[2] Ann. Phys. (Leipzig) 7 (1998), Proceedings of the Workshop on Superluminal(?) Velocities, ed. P. Mittelstaedt and G. 
Nimtz. 



14 



[3] The numerical value of u should not be overinterpreted. Only energy differences, and the differences between the corre- 
sponding frequencies are invariant with respect to shifts of the origin and physically significant. 
[4] K. W. H. Stevens, Eur. J. Phys. 1, 98 (1980); J. Phys. C: Solid State Phys. 16, 3649 (1983). 
[5] M. Biittiker and R. Landauer, Phys. Rev. Lett. 49, 1739 (1982). 
[6] M. Biittiker, Phys. Rev. B 27, 6178 (1983). 
[7] P. Moretti, Physica Scripta 45, 18 (1992) 

[8] A. Ranfangi, D. Mugnai, P. Fabeni and P. Pazzi, Physica Scripta 42, 508 (1990). 

[9] A. Ranfangi, D. Mugnai and A. Agresti, Phys. Lett. A 158, 161 (1991). 
[10] N. Teranishi, A. M. Kriman and D. K. Ferry, Superlattices and Microstructures, 3, 509 (1987). 
[11] A. P. Jauho nd M. Jonson, Superlattices and Microstructures 6, 303 (1989). 
[12] S. Brouard and J. G. Muga, Phys. Rev. A 54, 3055 (1996). 

[13] K. E. Oughstun and G. C. Sherman, Electromagnetic Pulse Propagatton in causal Dielectrics, (Springer, Berlin 1997). 
[14] M. Biittiker and H. Thomas, Ann. Phys. (Leipzig) 7, 602 (1998); Superlattices and Microstructures 23 781 (1998). 
[15] Handbook of Mathematical Functions, edited by M. Abramowitz and L A. Stegun (Dover, New York, 1972). 
[16] V. N. Faddeyeva and N. M. Terentev, Tables of the probability integral for complex argument (Pergamon Press, New York, 
1961). 

[17] N. Bleistein and R. Handelsman, Asymptotic Expansions of Integrals (Dover, New York, 1986). 
[18] Writting ip = |^|e"^ the local average instantaneous frequency is defined as |23] 

^ ' dt 

Note that in this work the convention of signs for Fourier transforms and for Lb[x) is different from the one in [psf. 
[19] G. Diener, Ann. Phys. (Leipzig) 7, 639 (1998). 
[20] T. E. Hartman, J. Appl. Phys. 33, 3427 (1962). 
[21] S. Brouard, R. Sala and J. G. Muga, Phys. Rev. A 49, 4312 (1994). 
[22] V. Delgado, J. G. Muga, Ann. Phys. (NY) 248, 122 (1996). 
[23] L. Cohen, Time- Frequency analysis (Prentice Hall, New Jersey, 1995). 

[24] L. Cohen, Annals of the New York Academy of Sciences, 808 97 (1997); L. Cohen, "Characterization of transients", in 

SPIE Proceedings, Vol. 3069, to appear. 
[25] J. G. Muga, V. Delgado, R. Sala and R. F. Snider, J. Chem. Phys. 104, 7015 (1996). 



15 



Figure Captions 



1. Decimal logarithm of the amplified density, A = \tjj{x, t)\^ exp(2«;oa;), versus time for a signal emitted by a source 
located at x = with frequency ojq = 0.5, and observed at x = 7, 14, 21 (solid, dotted, and dottcd-dashcd lines). 
Noticeable for x = 7 is the contribution of the saddle point solution tps (dashed line) in comparison with the 
exact solution ip. At x = 14, 21 the exact wave function tp and the saddle point solution ips are indistinguishable 
on the scale of the figure. 



2. Evolution of the average instantaneous frequency ui, saddle point frequency Us, and signal frequency cjq versus 
time (solid, dashed and dotted lines) for x = 7 and ujq = 0.5. The circle marks the time ttr at which the pole 
contribution and saddle point contribution are equal in magnitude. 



3. Decimal logarithm of the amplified density, A = \xjj{x,t)\^ cxp{2kqx), versus time for x = 0.05 and x = 0.45 
(solid and dashed thick lines respectively). The approximations provided by the saddle point solution tps are 
also shown (solid and dashed thin lines), = 0.5. The empty and filled squares mark the values of the traversal 
time T for a; = 0.05 and x = 0.45 respectively. 



4. Decimal logarithm of the amplified density, A = \tp{x, t)]"^ exp{2Kox), for the frequency band limited source and 
of the approximate solution ID^P versus time (solid and dashed lines respectively), for a; = 50 (lower set), and 
X = 135 (upper set), wq = 0.5, Aw = 0.12. 



5. Decimal logarithm of the amplified density, A = \ip{x,t)\'^ exp{2Kox), of the approximations [f+P and 

versus time (solid, dashed and dotted lines respectively), for x = 13.5. The square marks the value of r. 
wo = 0.5, Acu = 0.12 



6. Frequency evolution uj, Wg, and versus time (solid, dashed and dotted lines) of the frequency band limited 
source for x = 13.5, = 0.5, and Aw = 0.12 



7. Contour map of the spectrogram for wo = 0.5, T = 52.36, x = 135. Notice the transition from the kinetic to the 
potential regime and the side wings around the main peak. 



8. Spectrogram S{t, w = 0, 5) for wq = 0.5, T = 52.36, x = 135. The square marks the value of r = 95.4. (This is 
the period that corresponds to Aw = 0.12 in Figure D) Tq = 12.57. 



9. Spectrogram S{t,LU = 0,5) for wq = 0.5, T = 12.57, x = 135. The square marks the value of r = 95.4. 
To = 12.57. 
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